Bioinformatics and system biology approach to identify the influences among COVID-19, influenza, and HIV on the regulation of gene expression.

Background: Coronavirus disease (COVID-19), caused by SARS-CoV-2, has emerged as a infectious disease, coexisting with widespread seasonal and sporadic influenza epidemics globally. Individuals living with HIV, characterized by compromised immune systems, face an elevated risk of severe outcomes and increased mortality when affected by COVID-19. Despite this connection, the molecular intricacies linking COVID-19, influenza, and HIV remain unclear. Our research endeavors to elucidate the shared pathways and molecular markers in individuals with HIV concurrently infected with COVID-19 and influenza. Furthermore, we aim to identify potential medications that may prove beneficial in managing these three interconnected illnesses. Methods: Sequencing data for COVID-19 (GSE157103), influenza (GSE185576), and HIV (GSE195434) were retrieved from the GEO database. Commonly expressed differentially expressed genes (DEGs) were identified across the three datasets, followed by immune infiltration analysis and diagnostic ROC analysis on the DEGs. Functional enrichment analysis was performed using GO/KEGG and Gene Set Enrichment Analysis (GSEA). Hub genes were screened through a Protein-Protein Interaction networks (PPIs) analysis among DEGs. Analysis of miRNAs, transcription factors, drug chemicals, diseases, and RNA-binding proteins was conducted based on the identified hub genes. Finally, quantitative PCR (qPCR) expression verification was undertaken for selected hub genes. Results: The analysis of the three datasets revealed a total of 22 shared DEGs, with the majority exhibiting an area under the curve value exceeding 0.7. Functional enrichment analysis with GO/KEGG and GSEA primarily highlighted signaling pathways associated with ribosomes and tumors. The ten identified hub genes included IFI44L, IFI44, RSAD2, ISG15, IFIT3, OAS1, EIF2AK2, IFI27, OASL, and EPSTI1. Additionally, five crucial miRNAs (hsa-miR-8060, hsa-miR-6890-5p, hsa-miR-5003-3p, hsa-miR-6893-3p, and hsa-miR-6069), five essential transcription factors (CREB1, CEBPB, EGR1, EP300, and IRF1), and the top ten significant drug chemicals (estradiol, progesterone, tretinoin, calcitriol, fluorouracil, methotrexate, lipopolysaccharide, valproic acid, silicon dioxide, cyclosporine) were identified. Conclusion: This research provides valuable insights into shared molecular targets, signaling pathways, drug chemicals, and potential biomarkers for individuals facing the complex intersection of COVID-19, influenza, and HIV. These findings hold promise for enhancing the precision of diagnosis and treatment for individuals with HIV co-infected with COVID-19 and influenza.

Computational identification and experimental verification of a novel signature based on SARS-CoV-2-related genes for predicting prognosis, immune microenvironment and therapeutic strategies in lung adenocarcinoma patients.

Background: Early research indicates that cancer patients are more vulnerable to adverse outcomes and mortality when infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Nonetheless, the specific attributes of SARS-CoV-2 in lung Adenocarcinoma (LUAD) have not been extensively and methodically examined. Methods: We acquired 322 SARS-CoV-2 infection-related genes (CRGs) from the Human Protein Atlas database. Using an integrative machine learning approach with 10 algorithms, we developed a SARS-CoV-2 score (Cov-2S) signature across The Cancer Genome Atlas and datasets GSE72094, GSE68465, and GSE31210. Comprehensive multi-omics analysis, including assessments of genetic mutations and copy number variations, was conducted to deepen our understanding of the prognosis signature. We also analyzed the response of different Cov-2S subgroups to immunotherapy and identified targeted drugs for these subgroups, advancing personalized medicine strategies. The expression of Cov-2S genes was confirmed through qRT-PCR, with GGH emerging as a critical gene for further functional studies to elucidate its role in LUAD. Results: Out of 34 differentially expressed CRGs identified, 16 correlated with overall survival. We utilized 10 machine learning algorithms, creating 101 combinations, and selected the RFS as the optimal algorithm for constructing a Cov-2S based on the average C-index across four cohorts. This was achieved after integrating several essential clinicopathological features and 58 established signatures. We observed significant differences in biological functions and immune cell statuses within the tumor microenvironments of high and low Cov-2S groups. Notably, patients with a lower Cov-2S showed enhanced sensitivity to immunotherapy. We also identified five potential drugs targeting Cov-2S. In vitro experiments revealed a significant upregulation of GGH in LUAD, and its knockdown markedly inhibited tumor cell proliferation, migration, and invasion. Conclusion: Our research has pioneered the development of a consensus Cov-2S signature by employing an innovative approach with 10 machine learning algorithms for LUAD. Cov-2S reliably forecasts the prognosis, mirrors the tumor's local immune condition, and supports clinical decision-making in tumor therapies.

Probing molecular pathways: Illuminating the connection between COVID-19 and Alzheimer's disease through the endocannabinoid system dynamics.

Our study investigates the molecular link between COVID-19 and Alzheimer's disease (AD). We aim to elucidate the mechanisms by which COVID-19 may influence the onset or progression of AD. Using bioinformatic tools, we analyzed gene expression datasets from the Gene Expression Omnibus (GEO) database, including GSE147507, GSE12685, and GSE26927. Intersection analysis was utilized to identify common differentially expressed genes (CDEGs) and their shared biological pathways. Consensus clustering was conducted to group AD patients based on gene expression, followed by an analysis of the immune microenvironment and variations in shared pathway activities between clusters. Additionally, we identified transcription factor-binding sites shared by CDEGs and genes in the common pathway. The activity of the pathway and the expression levels of the CDEGs were validated using GSE164805 and GSE48350 datasets. Six CDEGs (MAL2, NECAB1, SH3GL2, EPB41L3, MEF2C, and NRGN) were identified, along with a downregulated pathway, the endocannabinoid (ECS) signaling pathway, common to both AD and COVID-19. These CDEGs showed a significant correlation with ECS activity (p < 0.05) and immune functions. The ECS pathway was enriched in healthy individuals' brains and downregulated in AD patients. Validation using GSE164805 and GSE48350 datasets confirmed the differential expression of these genes in COVID-19 and AD tissues. Our findings reveal a potential pathogenetic link between COVID-19 and AD, mediated by CDEGs and the ECS pathway. However, further research and multicenter evidence are needed to translate these findings into clinical applications.

Properties and Mechanisms of Deletions, Insertions, and Substitutions in the Evolutionary History of SARS-CoV-2.

SARS-CoV-2 has accumulated many mutations since its emergence in late 2019. Nucleotide substitutions leading to amino acid replacements constitute the primary material for natural selection. Insertions, deletions, and substitutions appear to be critical for coronavirus's macro- and microevolution. Understanding the molecular mechanisms of mutations in the mutational hotspots (positions, loci with recurrent mutations, and nucleotide context) is important for disentangling roles of mutagenesis and selection. In the SARS-CoV-2 genome, deletions and insertions are frequently associated with repetitive sequences, whereas C>U substitutions are often surrounded by nucleotides resembling the APOBEC mutable motifs. We describe various approaches to mutation spectra analyses, including the context features of RNAs that are likely to be involved in the generation of recurrent mutations. We also discuss the interplay between mutations and natural selection as a complex evolutionary trend. The substantial variability and complexity of pipelines for the reconstruction of mutations and the huge number of genomic sequences are major problems for the analyses of mutations in the SARS-CoV-2 genome. As a solution, we advocate for the development of a centralized database of predicted mutations, which needs to be updated on a regular basis.

Inflammasome-Related Genetic Polymorphisms as Severity Biomarkers of COVID-19.

The most critical forms of coronavirus disease 2019 (COVID-19) are associated with excessive activation of the inflammasome. Despite the COVID-19 impact on public health, we still do not fully understand the mechanisms by which the inflammatory response influences disease prognosis. Accordingly, we aimed to elucidate the role of polymorphisms in the key genes of the formation and signaling of the inflammasome as biomarkers of COVID-19 severity. For this purpose, a large and well-defined cohort of 377 COVID-19 patients with mild (n = 72), moderate (n = 84), severe (n = 100), and critical (n = 121) infections were included. A total of 24 polymorphisms located in inflammasome-related genes (NLRP3, NLRC4, NLRP1, CARD8, CASP1, IL1B, IL18, NFKB1, ATG16L1, and MIF) were genotyped in all of the patients and in the 192 healthy controls (HCs) (who were without COVID-19 at the time of and before the study) by RT-qPCR. Our results showed that patients with mild, moderate, severe, and critical COVID-19 presented similar allelic and genotypic distribution in all the variants studied. No statistically significant differences in the haplotypic distribution of NLRP3, NLRC4, NLRP1, CARD8, CASP1, IL1B, and ATG16L1 were observed between COVID-19 patients, who were stratified by disease severity. Each stratified group of patients presented a similar genetic distribution to the HCs. In conclusion, our results suggest that the inflammasome polymorphisms studied are not associated with the worsening of COVID-19.

Next-generation sequencing of host genetics risk factors associated with COVID-19 severity and long-COVID in Colombian population.

Coronavirus disease 2019 (COVID-19) was considered a major public health burden worldwide. Multiple studies have shown that susceptibility to severe infections and the development of long-term symptoms is significantly influenced by viral and host factors. These findings have highlighted the potential of host genetic markers to identify high-risk individuals and develop target interventions to reduce morbimortality. Despite its importance, genetic host factors remain largely understudied in Latin-American populations. Using a case-control design and a custom next-generation sequencing (NGS) panel encompassing 81 genetic variants and 74 genes previously associated with COVID-19 severity and long-COVID, we analyzed 56 individuals with asymptomatic or mild COVID-19 and 56 severe and critical cases. In agreement with previous studies, our results support the association between several clinical variables, including male sex, obesity and common symptoms like cough and dyspnea, and severe COVID-19. Remarkably, thirteen genetic variants showed an association with COVID-19 severity. Among these variants, rs11385942 (p < 0.01; OR = 10.88; 95% CI = 1.36-86.51) located in the LZTFL1 gene, and rs35775079 (p = 0.02; OR = 8.53; 95% CI = 1.05-69.45) located in CCR3 showed the strongest associations. Various respiratory and systemic symptoms, along with the rs8178521 variant (p < 0.01; OR = 2.51; 95% CI = 1.27-4.94) in the IL10RB gene, were significantly associated with the presence of long-COVID. The results of the predictive model comparison showed that the mixed model, which incorporates genetic and non-genetic variables, outperforms clinical and genetic models. To our knowledge, this is the first study in Colombia and Latin-America proposing a predictive model for COVID-19 severity and long-COVID based on genomic analysis. Our study highlights the usefulness of genomic approaches to studying host genetic risk factors in specific populations. The methodology used allowed us to validate several genetic variants previously associated with COVID-19 severity and long-COVID. Finally, the integrated model illustrates the importance of considering genetic factors in precision medicine of infectious diseases.

SARS-CoV-2 variant with the spike protein mutation F306L in the southern border provinces of Thailand.

The southernmost part of Thailand is a unique and culturally diverse region that has been greatly affected by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) outbreak during the coronavirus disease-2019 pandemic. To gain insights into this situation, we analyzed 1942 whole-genome sequences of SARS-CoV-2 obtained from the five southernmost provinces of Thailand between April 2021 and March 2022, together with those publicly available in the Global Initiative on Sharing All Influenza Data database. Our analysis revealed evidence for transboundary transmissions of the virus in and out of the five southernmost provinces during the study period, from both domestic and international sources. The most prevalent viral variant in our sequence dataset was the Delta B.1.617.2.85 variant, also known as the Delta AY.85 variant, with many samples carrying a non-synonymous mutation F306L in their spike protein. Protein-protein docking and binding interface analyses suggested that the mutation may enhance the binding between the spike protein and host cell receptor protein angiotensin-converting enzyme 2, and we found that the mutation was significantly associated with an increased fatality rate. This mutation has also been observed in other SARS-CoV-2 variants, suggesting that it is of particular interest and should be monitored.

The fatal contribution of serine protease-related genetic variants to COVID-19 outcomes.

Introduction: Serine proteases play a critical role during SARS-CoV-2 infection. Therefore, polymorphisms of transmembrane protease serine 2 (TMPRSS2) and serpine family E member 1 (SERPINE1) could help to elucidate the contribution of variability to COVID-19 outcomes. Methods: To evaluate the genetic variants of the genes previously associated with COVID-19 outcomes, we performed a cross-sectional study in which 1536 SARS-CoV-2-positive participants were enrolled. TMPRSS2 (rs2070788, rs75603675, rs12329760) and SERPINE1 (rs2227631, rs2227667, rs2070682, rs2227692) were genotyped using the Open Array Platform. The association of polymorphisms with disease outcomes was determined by logistic regression analysis adjusted for covariates (age, sex, hypertension, type 2 diabetes, and obesity). Results: According to our codominant model, the GA genotype of rs2227667 (OR=0.55; 95% CI = 0.36-0.84; p=0.006) and the AG genotype of rs2227667 (OR=0.59; 95% CI = 0.38-0.91; p=0.02) of SERPINE1 played a protective role against disease. However, the rs2227692 T allele and TT genotype SERPINE1 (OR=1.45; 95% CI = 1.11-1.91; p=0.006; OR=2.08; 95% CI = 1.22-3.57; p=0.007; respectively) were associated with a decreased risk of death. Similarly, the rs75603675 AA genotype TMPRSS2 had an OR of 1.97 (95% CI = 1.07-3.6; p=0.03) for deceased patients. Finally, the rs2227692 T allele SERPINE1 was associated with increased D-dimer levels (OR=1.24; 95% CI = 1.03-1.48; p=0.02). Discussion: Our data suggest that the rs75603675 TMPRSS2 and rs2227692 SERPINE1 polymorphisms are associated with a poor outcome. Additionally, rs2227692 SERPINE1 could participate in hypercoagulable conditions in critical COVID-19 patients, and this genetic variant could contribute to the identification of new pharmacological targets and treatment strategies to block the inhibition of TMPRSS2 entry into SARS-CoV-2.

Comprehensive analysis of immunogenic cell death-related gene and construction of prediction model based on WGCNA and multiple machine learning in severe COVID-19.

The death of coronavirus disease 2019 (COVID-19) is primarily due to from critically ill patients, especially from ARDS complications caused by SARS-CoV-2. Therefore, it is essential to contribute an in-depth understanding of the pathogenesis of the disease and to identify biomarkers for predicting critically ill patients at the molecular level. Immunogenic cell death (ICD), as a specific variant of regulatory cell death driven by stress, can induce adaptive immune responses against cell death antigens in the host. Studies have confirmed that both innate and adaptive immune pathways are involved in the pathogenesis of SARS-CoV-2 infection. However, the role of ICD in the pathogenesis of severe COVID-19 has rarely been explored. In this study, we systematically evaluated the role of ICD-related genes in COVID-19. We conducted consensus clustering, immune infiltration analysis, and functional enrichment analysis based on ICD differentially expressed genes. The results showed that immune infiltration characteristics were altered in severe and non-severe COVID-19. In addition, we used multiple machine learning methods to screen for five risk genes (KLF5, NSUN7, APH1B, GRB10 and CD4), which are used to predict COVID-19 severity. Finally, we constructed a nomogram to predict the risk of severe COVID-19 based on the classification and recognition model, and validated the model with external data sets. This study provides a valuable direction for the exploration of the pathogenesis and progress of COVID-19, and helps in the early identification of severe cases of COVID-19 to reduce mortality.

Novel genetic association of the Furin gene polymorphism rs1981458 with COVID-19 severity among Indian populations.

SARS CoV-2, the causative agent for the ongoing COVID-19 pandemic, it enters the host cell by activating the ACE2 receptor with the help of two proteasesi.e., Furin and TMPRSS2. Therefore, variations in these genes may account for differential susceptibility and severity between populations. Previous studies have shown that the role of ACE2 and TMPRSS2 gene variants in understanding COVID-19 susceptibility among Indian populations. Nevertheless, a knowledge gap exists concerning the COVID-19 susceptibility of Furin gene variants among diverse South Asian ethnic groups. Investigating the role of Furin gene variants and their global phylogeographic structure is essential to comprehensively understanding COVID-19 susceptibility in these populations. We have used 450 samples from diverse Indian states and performed linear regression to analyse the Furin gene variant's with COVID-19 Case Fatality Rate (CFR) that could be epidemiologically associated with disease severity outcomes. Associated genetic variants were further evaluated for their expression and regulatory potential through various Insilco analyses. Additionally, we examined the Furin gene using next-generation sequencing (NGS) data from 393 diverse global samples, with a particular emphasis on South Asia, to investigate its Phylogeographic structure among diverse world populations. We found a significant positive association for the SNP rs1981458 with COVID-19 CFR (p < 0.05) among diverse Indian populations at different timelines of the first and second waves. Further, QTL and other regulatory analyses showed various significant associations for positive regulatory roles of rs1981458 and Furin gene, mainly in Immune cells and virus infection process, highlighting their role in host immunity and viral assembly and processing. The Furin protein-protein interaction suggested that COVID-19 may contribute to Pulmonary arterial hypertension via a typical inflammation mechanism. The phylogeographic architecture of the Furin gene demonstrated a closer genetic affinity of South Asia with West Eurasian populations. Therefore, it is worth proposing that for the Furin gene, the COVID-19 susceptibility of South Asians will be more similar to the West Eurasian population. Our previous studies on the ACE2 and TMPRSS2 genes showed genetic affinity of South Asian with East Eurasians and West Eurasians, respectively. Therefore, with the collective information from these three important genes (ACE2, TMPRSS2 and Furin) we modelled COVID-19 susceptibilityof South Asia in between these two major ancestries with an inclination towards West Eurasia. In conclusion, this study, for the first time, concluded the role of rs1981458 in COVID-19 severity among the Indian population and outlined its regulatory potential.This study also highlights that the genetic structure for COVID-19 susceptibilityof South Asia is distinct, however, inclined to the West Eurasian population. We believe this insight may be utilised as a genetic biomarker to identify vulnerable populations, which might be directly relevant for developing policies and allocating resources more effectively during an epidemic.

Unveiling hidden structural patterns in the SARS-CoV-2 genome: Computational insights and comparative analysis.

SARS-CoV-2, the causative agent of COVID-19, is known to exhibit secondary structures in its 5' and 3' untranslated regions, along with the frameshifting stimulatory element situated between ORF1a and 1b. To identify additional regions containing conserved structures, we utilized a multiple sequence alignment with related coronaviruses as a starting point. We applied a computational pipeline developed for identifying non-coding RNA elements. Our pipeline employed three different RNA structural prediction approaches. We identified forty genomic regions likely to harbor structures, with ten of them showing three-way consensus substructure predictions among our predictive utilities. We conducted intracomparisons of the predictive utilities within the pipeline and intercomparisons with four previously published SARS-CoV-2 structural datasets. While there was limited agreement on the precise structure, different approaches seemed to converge on regions likely to contain structures in the viral genome. By comparing and combining various computational approaches, we can predict regions most likely to form structures, as well as a probable structure or ensemble of structures. These predictions can be used to guide surveillance, prophylactic measures, or therapeutic efforts. Data and scripts employed in this study may be found at https://doi.org/10.5281/zenodo.8298680.

Reassessing human MHC-I genetic diversity in T cell studies.

The Major Histocompatibility Complex class I (MHC-I) system plays a vital role in immune responses by presenting antigens to T cells. Allele specific technologies, including recombinant MHC-I technologies, have been extensively used in T cell analyses for COVID-19 patients and are currently used in the development of immunotherapies for cancer. However, the immense diversity of MHC-I alleles presents challenges. The genetic diversity serves as the foundation of personalized medicine, yet it also poses a potential risk of exacerbating healthcare disparities based on MHC-I alleles. To assess potential biases, we analysed (pre)clinical publications focusing on COVID-19 studies and T cell receptor (TCR)-based clinical trials. Our findings reveal an underrepresentation of MHC-I alleles associated with Asian, Australian, and African descent. Ensuring diverse representation is vital for advancing personalized medicine and global healthcare equity, transcending genetic diversity. Addressing this disparity is essential to unlock the full potential of T cells for enhancing diagnosis and treatment across all individuals.

Metacell-based differential expression analysis identifies cell type specific temporal gene response programs in COVID-19 patient PBMCs.

By profiling gene expression in individual cells, single-cell RNA-sequencing (scRNA-seq) can resolve cellular heterogeneity and cell-type gene expression dynamics. Its application to time-series samples can identify temporal gene programs active in different cell types, for example, immune cells' responses to viral infection. However, current scRNA-seq analysis has limitations. One is the low number of genes detected per cell. The second is insufficient replicates (often 1-2) due to high experimental cost. The third lies in the data analysis-treating individual cells as independent measurements leads to inflated statistics. To address these, we explore a new computational framework, specifically whether "metacells" constructed to maintain cellular heterogeneity within individual cell types (or clusters) can be used as "replicates" for increasing statistical rigor. Toward this, we applied SEACells to a time-series scRNA-seq dataset from peripheral blood mononuclear cells (PBMCs) after SARS-CoV-2 infection to construct metacells, and used them in maSigPro for quadratic regression to find significantly differentially expressed genes (DEGs) over time, followed by clustering expression velocity trends. We showed that such metacells retained greater expression variances and produced more biologically meaningful DEGs compared to either metacells generated randomly or from simple pseudobulk methods. More specifically, this approach correctly identified the known ISG15 interferon response program in almost all PBMC cell types and many DEGs enriched in the previously defined SARS-CoV-2 infection response pathway. It also uncovered additional and more cell type-specific temporal gene expression programs. Overall, our results demonstrate that the metacell-pseudoreplicate strategy could potentially overcome the limitation of 1-2 replicates.

Polymorphisms of IFN signaling genes and FOXP4 influence the severity of COVID-19.

BACKGROUND: The clinical manifestations of COVID-19 range from asymptomatic, mild to moderate, severe, and critical disease. Host genetic variants were recognized to affect the disease severity. However, the genetic landscape differs among various populations. Therefore, we explored the variants associated with COVID-19 severity in the Guangdong population. METHODS: A total of 314 subjects were selected, of which the severe and critical COVID-19 patients were defined as "cases", and the mild and moderate patients were defined as "control". Twenty-two variants in interferon-related genes and FOXP4 were genotyped using the MassARRAY technology platform. RESULTS: IFN signaling gene MX1 rs17000900 CA + AA genotype was correlated with a reduced risk of severe COVID-19 in males (P = 0.001, OR = 0.050, 95%CI = 0.008-0.316). The AT haplotype comprised of MX1 rs17000900 and rs2071430 was more likely to protect against COVID-19 severity (P = 6.3E-03). FOXP4 rs1886814 CC genotype (P = 0.001, OR = 3.747, 95%CI = 1.746-8.043) and rs2894439 GA + AA genotype (P = 0.001, OR = 5.703, 95% CI = 2.045-15.903) were correlated with increased risk of severe COVID-19. Haplotype CA comprised of rs1886814 and rs2894439 was found to be correlated with adverse outcomes (P = 7.0E-04). FOXP4 rs1886814 CC (P = 0.0004) and rs2894439 GA + AA carriers had higher neutralizing antibody titers (P = 0.0018). The CA + AA genotype of MX1 rs17000900 tended to be correlated with lower neutralizing antibody titers than CC genotype (P = 0.0663), but the difference was not statistically significant. CONCLUSION: Our study found a possible association between MX1 and FOXP4 polymorphisms and the severity of COVID-19. Distinguishing high-risk patients who develop severe COVID-19 will provide clues for early intervention and individual treatment strategies.

The association between COVID-19 and infertility: Mendelian randomization analysis.

Since December 2019, COVID-19 has triggered a global pandemic. The association of COVID-19 with the long-term reproductive situation of women and males is not clear. Thus, our aim was to assess the causal association between COVID-19 and infertility using Mendelian randomization (MR) analysis based on the OpenGWAS database. Two-sample MR analysis was conducted using one genome-wide association study (GWAS) on COVID-19 and infertility in individuals of European ancestry. The summary data of genetic variation come from the GWAS in European populations. We applied several MR methods, including MR Egger, weighted median, inverse variance weighted, simple mode, weighted mode, to test causal relationships. After observing the statistical analysis results of MR, we conducted sensitivity analysis to test robustness. After gene prediction, it was found that there was no clear causal relationship between COVID-19 and male infertility in MR analysis [OR 0.4702 (95% CI, 0.1569-1.4093), P = .178]. Moreover, COVID-19 was not associated with female infertility [OR 0.9981 (95% CI, 0.763-1.544), P = .646]. Sensitivity analysis showed that the MR results were robust [level pleiotropy, male: (MR-Egger, intercept = 0.1967434; se = 0.1186876; P = .2392406); female: (MR-Egger, intercept = -0.05902506; se = 0.05362049; P = .3211367)]. To further validate the impact of COVID-19 on infertility, we added a covariate (sex hormone binding global levels, abortion) to the MR analysis, which is a multivariate MR analysis. According to univariate and multivariate MR analyses, the evidence does not support that COVID-19 is a causal risk factor for infertility in European population. This information can provide information for doctors in reproductive centers when managing infertility patients.

Molecular Profile of Variants Potentially Associated with Severe Forms of COVID-19 in Amazonian Indigenous Populations.

Coronavirus disease 2019 (COVID-19) is an infection caused by SARS-CoV-2. Genome-wide association studies (GWASs) have suggested a strong association of genetic factors with the severity of the disease. However, many of these studies have been completed in European populations, and little is known about the genetic variability of indigenous peoples' underlying infection by SARS-CoV-2. The objective of the study is to investigate genetic variants present in the genes AQP3, ARHGAP27, ELF5L, IFNAR2, LIMD1, OAS1 and UPK1A, selected due to their association with the severity of COVID-19, in a sample of indigenous people from the Brazilian Amazon in order to describe potential new and already studied variants. We performed the complete sequencing of the exome of 64 healthy indigenous people from the Brazilian Amazon. The allele frequency data of the population were compared with data from other continental populations. A total of 66 variants present in the seven genes studied were identified, including a variant with a high impact on the ARHGAP27 gene (rs201721078) and three new variants located in the Amazon Indigenous populations (INDG) present in the AQP3, IFNAR2 and LIMD1 genes, with low, moderate and modifier impact, respectively.

Evolution of nasal and olfactory infection characteristics of SARS-CoV-2 variants.

SARS-CoV-2 infection of the upper airway and the subsequent immune response are early, critical factors in COVID-19 pathogenesis. By studying infection of human biopsies in vitro and in a hamster model in vivo, we demonstrated a transition in nasal tropism from olfactory to respiratory epithelium as the virus evolved. Analyzing each variant revealed that SARS-CoV-2 WA1 or Delta infect a proportion of olfactory neurons in addition to the primary target sustentacular cells. The Delta variant possessed broader cellular invasion capacity into the submucosa, while Omicron displayed enhanced nasal respiratory infection and longer retention in the sinonasal epithelium. The olfactory neuronal infection by WA1 and the subsequent olfactory bulb transport via axon were more pronounced in younger hosts. In addition, the observed viral clearance delay and phagocytic dysfunction in aged olfactory mucosa were accompanied by a decline of phagocytosis-related genes. Further, robust basal stem cell activation contributed to neuroepithelial regeneration and restored ACE2 expression postinfection. Together, our study characterized the nasal tropism of SARS-CoV-2 strains, immune clearance, and regeneration after infection. The shifting characteristics of viral infection at the airway portal provide insight into the variability of COVID-19 clinical features, particularly long COVID, and may suggest differing strategies for early local intervention.

PANoptosis, an indicator of COVID-19 severity and outcomes.

Coronavirus disease 2019 (COVID-19) has been wreaking havoc for 3 years. PANoptosis, a distinct and physiologically relevant inflammatory programmed cell death, perpetuates cytokine storm and multi-organ injuries in COVID-19. Although PANoptosis performs indispensable roles in host defense, further investigation is needed to elucidate the exact processes through which PANoptosis modulates immunological responses and prognosis in COVID-19. This study conducted a bioinformatics analysis of online single-cell RNA sequence (scRNA-seq) and bulk RNA-seq datasets to explore the potential of PANoptosis as an indicator of COVID-19 severity. The degree of PANoptosis in bronchoalveolar lavage fluid (BALF) and peripheral blood mononuclear cells (PBMC) indicated the severity of COVID-19. Single-cell transcriptomics identified pro-inflammatory monocytes as one of the primary sites of PANoptosis in COVID-19. The study subsequently demonstrated the immune and metabolic characteristics of this group of pro-inflammatory monocytes. In addition, the analysis illustrated that dexamethasone was likely to alleviate inflammation in COVID-19 by mitigating PANoptosis. Finally, the study showed that the PANoptosis-related genes could predict the intensive care unit admission (ICU) and outcomes of COVID-19 patients who are hospitalized.

Bulk and Single-Cell RNA Sequencing Elucidate the Etiology of Severe COVID-19.

Coronavirus disease 2019 (COVID-19) is a type of pneumonia caused by the SARS-CoV-2 coronavirus. It can cause acute pulmonary and systemic inflammation, which can lead to death in severely ill patients. This study explores the potential reasons behind severe COVID-19 and its similarities to systemic autoimmune diseases. This study reviewed unbiased high-throughput gene expression datasets, including next-generation and single-cell RNA sequencing. A total of 27 studies and eight meta-analyses were reviewed. The studies indicated that severe COVID-19 is associated with the upregulation of genes involved in pro-inflammatory, interferon, and cytokine/chemokine pathways. Additionally, changes were observed in the proportions of immune cell types in the blood and tissues, along with degenerative alterations in lung epithelial cells. Genomic evidence also supports the association of severe COVID-19 with various inflammatory syndromes, such as neuronal COVID-19, acute respiratory distress syndrome, vascular inflammation, and multisystem inflammatory syndrome. In conclusion, this study suggests that gene expression profiling plays a significant role in elucidating the etiology of severe COVID-19.